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Abstract. We analyze the large-distance behaviour of static quark-anti-quark pair 
■ correlations in QCD. The singlet free energy is calculated and the entropy contribution 

to it is identified allowing us to calculate the excess internal energy. The free energy 
has a sharp drop in the critical region, leading to sharp peaks in both excess entropy 
and internal energy. 
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1. Static quark anti-quark free energies 



Among its many other applications, finite temperature lattice QCD is a useful tool to 
study the modification of quark-anti-quark forces n a hot medium. Recent progress in 
large-scale simulations allowed us to address this problem in QCD with realistic quark 
masses. This study is based on a ongoing RBC-Bielefeld joint research effort which is 
fully described in Ref.[2]. Here we will note that calculations are performed at physical 
Q-i! strange quark mass and two light quark masses m„ = 0.1m s and 0.2m s , with m s being 
the strange quark mass; two cutoffs, corresponding to temporal extent N t = 4, 6 were 
used. 

We discuss here the asymptotic large distance behavior of qq renormalized free 
energies. This is a gauge invariant observable. However, it is most suitable to extract 
it from a qq singlet free energy, which is gauge variant and thus has been calculated by 
us in Coulomb gauge [SI El El El El CEDl ED] • To extract singlet contribution to the free 
energy one uses a correlation function of two static quarks, represented by temporal 
Wilson lines, one representing a static quark at x = and another - a static anti-quark 
at distance r, 

exp(— Fx(r, T)/T + C) = -(TrW(r)W j (0)) . (1) 
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This is, of course, not the full free energy of the system but the excess free energy 
which is generated by adding a gg-pair to the hot bath. In the zero temperature limit 
the singlet free energy defined above coincides with the well known static potential, as 
demonstrated in e.g. Ref. [12J. At finite temperature Fi(r,T) gives information about 
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in-medium modification of inter-quark forces and color screening. Both the singlet free 
energy and the zero temperature static potential are defined up to an additive constant 
C which depends on the lattice cutoff. Our renormalization scheme consists of using 
the zero temperature potential as normalization at short distances, where the effect of 
the medium is negligible. In the transition region we used a Cornell-type ansatz for the 
potential, 

0.385 1.263 

V{r) = -— + -3-r • (2) 
' 'o 

which gives a good parametrization of the lattice data for the zero temperature static 
potential calculated in Ref . [2] . At much higher temperatures, however, as we probe the 
free energy at increasingly shorter distances, where the effects of the running coupling 
can become important, we calculated V(r) explicitly at every temperature and used 
them for normalization. Our results for the singlet free energy and the scaling study 
have been presented in Ref. p. Here we will focus on the large-distance quantities. 



2. Renormalized Polyakov Loop 

The expectation value of the Polyakov loop, the trace of temporal Wilson line, (L(f)) = 
(TrW(r)), is the order parameter for the deconfining transition in pure gauge theories. 
In QCD dynamical quarks the relevant Z(3) symmetry is broken explicitly and (L(f)) 
thus is no longer an order parameter. Still it remains an interesting quantity to study the 
deconfinement transition as it shows a rapid increase in the crossover region [21 [13j [TJ] 
and is one of the ways which can be used to determine the transition temperature [21 IT3]. 
The unrenormalized Polyakov loop defined above strongly depends on the lattice cutoff 
and has no meaningful continuum limit. On the other hand a correlator of Polyakov 
loops is a physical quantity and corresponds to the color averaged free energy up to a 
normalization constant. Due to the cluster decomposition of correlation functions we 
have, 

exp(-F(r,T) + C) = l(L(filS(p))\™ = \(L(0))\ 2 - (3) 

Here the normalization constant C is the same as in EqJTJ Moreover, at large distances 
the relative color orientation of the gg-pair is not important and thus the color singlet 
free energy and the color averaged free energy approach the same constant value F^iT). 
Therefore, following Ref. [6], we define the renormalized Polyakov loop as 

L ren (T)=exp(-^^). (4) 

First, in order to check the scaling behaviour, let us concentrate on the transition region. 
As one can see from FigJT] L ren (T) shows an almost universal behavior as function of 
T/T c for all quark masses studied by us, including the case of three degenerate flavors. 
This suggests that in the region of the small quark masses studied here, the flavor and 
quark mass dependence of the free energy can almost entirely be absorbed in the flavor 
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Figure 1. The renormalized Polyakov in the vicinity of the transition calculated on 
16 3 x 6, 16 3 x 4 lattices for different quark masses (left) and in the entire temperature 
range (right) 



and quark mass dependence of the transition temperature T c . Note that the results 
obtained on 16 3 x 4 lattice are in remarkably good agreement with the results obtained 
on 16 3 x 6 lattices, indicating again that cutoff effects are small. For comparison with 
results obtained using different simulation techniques we refer to Ref.p]. For higher 
temperatures we simulated a large ensemble of 16 3 x 4 (finite temperature). Here results 
obtained for the static potential at T = on f 6 3 x 32 (zero temperature) lattices where 
used for the normalization. Results for the renormalized Polyakov loop are presented 
on Fig{T] (right). As one can see it shows a rapid change in the temperature region from 
approximately f 80MeV to 230MeV and almost reaches a plateau at about 400MeV. 

3. Entropy and Internal Energy 

The excess free energy obviously has a behaviour similar to that of the renormalized 
Polyakov loop through which it was defined. Its variation with temperature defines the 
excess energy and entropy, respectively, derived from standard thermodynamic relations. 
As one can see from FigJS] the entropy and internal energy have sharp peaks in the 
transition region. This can be intuitively understood in a following way: while QCD 
with non-zero quark masses does not undergo phase transition the crossover is rapid 
enough to share some important properties, such as enhancement of the long distance 
correlations. Therefore, adding a gg-pair to such a system would require very large, 
possibly overlapping, gluonic clouds to screen the color charges, which, in turn, gives 
rise to the aforementioned peak in the internal energy. 

As a phenomenological consequences of this calculation we note that potential 
models which make use of our data(Refs.[16],[l7j) fail to reproduce more rigorous results 
obtained in Ref.[15j. namely charmonia survival till about 2.5T C with simultaneous 
dissolution of the \ c . This shows the need for more refined phenomenological approaches. 
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Figure 2. Excess entropy and internal energy of the quark-anti-quark pair calculated 
on 16 3 x 4 lattices for m = 0.1m s 



4. Conclusions 

We have calculated the singlet free energy of a static gg-pair in QCD with realistic quark 
masses. We have found that this and related quantities show little cut-off dependence 
and can be calculated reliably on relatively coarse lattices. The renormalized Polyakov 
loop derived from the large distance limit of the singlet free energy has been calculated 
and shown to be in the scaling regime. The free energy experiences rapid drop in the 
transition region, which leads to sharp peaks in the excess entropy and internal energy. 
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